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At finite temperature, chiral quark models do not incorporate large gauge invari- 
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ance which implies genuinely non-perturbative finite temperature gluonic degrees 
of freedom. Motivated by this observation, we describe how the coupling of the 
Polyakov loop as an independent degree of freedom to quarks not only accounts 
for large gauge invariance, but also allows to establish in a dynamical way the 
interaction between composite hadronic states such as Goldstone bosons to finite 
Q H ' temperature non-perturbative gluons in a medium which can undergo a confinement- 

deconfinement phase transition. 

1 Large Gauge Transformations 

X ' 

One feature of gauge theories like QCD at finite temperatures in the imaginary time 
formulation [1, 2, 3] is the non-perturbative manifestation of the non Abelian gauge sym- 
metry. In the Polyakov gauge, where d^A^ = and A4 is a diagonal and traceless N c x N c 
matrix, and N c is the number of colors, there is still some freedom in choosing the gluon 
field. Let us consider for instance the periodic gauge transformation [4, 5] 

g ( X4 ) = e *"**V/J , (1) 

where A is a color traceless diagonal matrix of integers. We call it a large gauge trans- 
formation (LGT) since it cannot be considered to be close to the identity 1 . The gauge 
transformation on the A4 component of the gluon field is 

At — > At + —A . (2) 
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1 Note that they are not large in the topological sense, as discussed in [4, 5]. 
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Thus, invariance under the LGT, Eq. (1), implies a constant shift in the A\ gluon ampli- 
tudes, meaning that A 4 is not uniquely defined by the Polyakov gauge condition. These 
ambiguities on the choice of the gauge field within a given gauge fixing are usually called 
Gribov copies. The requirement of gauge invariance actually implies identifying all ampli- 
tudes differing by a multiple of 2ir//3, which means periodicity in the diagonal amplitudes 
of A A of period 2tt/(3. Perturbation theory, which corresponds to expanding in powers of 
small A A fields manifestly breaks gauge invariance at finite temperature, because a Taylor 
expansion on a periodic function violates the periodicity behavior. Thus, taking into ac- 
count these Gribov replicas is equivalent to explicitly deal with genuine non-perturbative 
finite temperature gluonic degrees of freedom. A way of automatically taking into ac- 
count LGT is by considering the Polyakov loop Q as an independent variable, which in 
the Polyakov gauge becomes a diagonal unitary matrix 

n = e i/3A4(£) (3) 

invariant under the set of transformations given by Eq. (1). The relevance of the Polyakov 
loop in practical calculations is well recognized [1] but seldomly taken into account in 
high temperature calculations where large gauge invariance is manifestly broken since the 
gluon field is considered to be small. We have recently developed an expansion keeping 
these symmetries in general theories and applied it to QCD at the one quark+gluon loop 
level [6, 7]. 



2 The Center Symmetry 

In pure gluodynamics, or in the quenched approximation (valid for heavy quarks) at 
finite temperature there is actually a larger symmetry since one can extend the periodic 
transformations to aperiodic ones [3], 

g(x A + P)=zg(x A ), z N < = l (4) 

so that z is an element of the center Z(N C ) of the group SU(N C ). This center symmetry 
is a symmetry of the action as well as the gluon field boundary conditions. An example 
of such a transformation in the Polyakov gauge is given by 

On the A4 component of the gluon field produces 

A,^A, + —A. (6) 

Thus, in the quenched approximation the period is N c times smaller than in full QCD. 
Under these transformations the gluonic action, measure and boundary conditions are 
invariant. The Polyakov loop, however, transforms as the fundamental representation of 
the Z(N C ) group, i.e. f2 — > zQ, yielding (Q) = z{Q) and hence (Q) = in the unbroken 
center symmetry phase. At high temperatures one expects perturbation theory to hold, 
the gluon field amplitude becomes small and hence (Q) — > 1, justifying the choice of fl as 
an order parameter for a confinement-deconfinement phase transition. More generally, in 
the confining phase 

(Q n ) =0 for n ^ mN c (7) 
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with m an arbitrary integer. The antiperiodic quark fields at the end of the Euclidean 
imaginary interval transform as q(x, j3) = —q(x,0) — > zq(x,j3) = —q(x,0), so that the 
center symmetry is explicitly broken by the presence of dynamical quarks. A direct con- 
sequence of such a property is that, in the quenched approximation non-local condensates 
fulfill a selection rule of the form, 

(q(n(3)q(0)) = for n ^ mN c (8) 

since under the large aperiodic transformations given by Eq. (5) we have q{n[3)q(0) — > 
z~ n q(n/3)q(0). This selection rule has some impact on chiral quark models. 



3 Chiral quark models at finite temperature 

To fully appreciate the role played by the center symmetry in chiral quark models (for a 
recent review on such models see e.g. Ref. [8] and references therein) let us evaluate the 
chiral condensate at finite temperature. At the one loop level one has 2 

r rl 3 k 1 

(*>• = 4M^E/^ F ^ TCT (») 

where u) n = 2nT{n + 1/2) are the fermionic Matsubara frequencies, M is the constituent 
quark mass and Tr c stands for the color trace in the fundamental representation which in 
this case trivially yields a N c factor. Possible finite cut-off corrections, appearing in the 
chiral quark models such as the NJL model at finite temperature have been neglected. This 
is a reasonable approximation as long as the temperature is low enough T < A ~ lGeV. 
The condensate can be rewritten as 

<ot>* = E(-i) n (^Mo)> (io) 

n 

in terms of nonlocal Euclidean condensates at zero temperature. After Poisson resumma- 
tion, at low temperatures we have 

oo TM 2 
(qq)* = (qq)+%N c Y,(-) n —K 1 {Mn/T), 

n=l 71 



(,)-f(-)^m 3/2 e^ 
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{^^) e > ( U ) 



where the asymptotic expansion of the modified Bessel function K\ has been used. One 
can interpret the previous formula for the condensate in terms of statistical Boltzmann 
factors, since at large Euclidean coordinates the fermion propagator behaves as S(if3, x) ~ 
e~ M/3 , so that we have contributions from multiquark states. This is a problem since 
it means that the heat bath is made out of free constituent quarks without any color 
clustering 3 . Another problem comes from comparison with Chiral Perturbation Theory 



2 We use an asterisk to denote finite temperature observables. 

3 One could think that this is a natural consequence of the lack of confinement in chiral quark models 
such as NJL. Contrary to naive expectations this is not necessarily the case; Boltzmann factors occur in 
quark models with analytic confinement such as the Spectral Quark Model [9]. There the condensate is 
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at Finite Temperature [10]. In the chiral limit, i.e., for m n <C 2ttT <C 4:irf n the leading 
thermal corrections to the quark condensate are given by 

<??>* ChpT = (qq) \1 ~ - 3^ + ...J. (13) 

This formula is derived under the assumption that there is no temperature dependence 
of the low energy constants, i.e. L* ~ Lj so that the whole effect is due to thermal pion 
loops. Thus, the finite temperature correction is iV c -suppressed as compared to the zero 
temperature value. This is not what one sees in chiral quark model calculations; in the 
large N c limit there is a finite temperature correction, which would mean that the low 
energy constants which appear in the chiral Lagrangian would have a genuine tree level 
temperature dependence, L* — Li ~ N c e~ M / T . To obtain the ChPT result of Eq. (13) 
pion loops have to be considered [11] and dominate for T <C M. The problem is that 
already without pion loops chiral quark models predict a chiral phase transition at about 
T c ~ 170 MeV, in remarkable but perhaps unjustified agreement with lattice calculations. 



4 Coupling the Polyakov loop 

In the Polyakov gauge one can formally keep track of large gauge invariance at finite 
temperature by coupling gluons to the model in a minimal way. This means in practice 
using the modified fermionic Matsubara frequencies [4, 5] 

Co n = 2vrT(n + l/2 + z/), i/ = (27ri) _1 logft (14) 

which are shifted by the logarithm of the Polyakov loop which we assume for simplicity 
to be x independent. Previous work have coupled similarly fl on pure phenomenological 
grounds [12, 13, 14], but the key role played by the implementation of large gauge in- 
variance was not recognized. This is the only place where explicit dependence on colour 
degrees of freedom appear. This coupling introduces a colour source into the problem for 
a fixed A field and projection onto the colour neutral states by integrating over the A 
field, in a gauge invariant manner, as required. Actually, at the one quark loop level there 
is an accidental Z(N C ) symmetry in the model which generates a similar selection rule as 
in pure gluodynamics, from which a strong thermal suppression, 0(e~ NcM / T ) follows. In 
this way compliance with ChPT can be achieved since now L* — Lj ~ e - N c M / T but also 
puts some doubts on whether chiral quark models still predict a chiral phase transition 
at realistic temperatures. This question has been addressed using specific potentials for 
the Polyakov loop either based on one loop perturbation theory for massive gluons [13] in 
the high temperature approximation or on strong coupling expansions on the lattice [14] . 
In both cases similar mean field qualitative features are displayed; the low temperature 
evolution is extremely flat, but there appears a rapid change in the critical region, so that 
(qq)* ~ (qq) when (Q) ~ and (qq) ~ when (Q) ~ 1. A more general discussion and 

given by 

= tanh(Af/2T) = 1 - 2e- M / T + 2e- 2M / T + . . . (12) 
where M — Ms/2, despite the absence of poles in the quark propagator. 
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diagramatic interpretation of these issues as well as the influence of higher quark loop 
effects and dynamical Polyakov loop contributions will be presented elsewhere [15] pro- 
viding a justification of the one quark loop approximation at least at low temperatures. 
There one obtains that the Polyakov loop effect can be factored out as follows 4 

<??>* = E^ Tr c ((-^r)(^/%(o)>- (is) 

This result is consistent with applying the center symmetry selection rule, Eq. (8), to 
the Z(N C ) breaking condensate, Eq. (10), of the chiral quark model without Polyakov 
loops. If one now takes a suitable average on Polyakov loop configurations consistent 
with center symmetry, i.e., including for each such configuration all its Gribov replicas, 
Eq. (7) applies. Schematically, this yields 

<??>* ~ £<?>WMo))~£e-^ M / T (i6) 

n n 

in the confining phase. (In the above sums each term carries a weight coming from the 
Polyakov loop average and phase space factors.) On the other hand in the unconfining 
phase, where the center symmetry is spontaneously broken, the Polyakov loop is nearly 
unity and one recovers the standard chiral quark models results, without Polyakov loop 
coupling. 



5 Chiral Lagrangians at finite temperature 

It is interesting to construct the coupling of Polyakov loops with composite pion fields 
at finite temperature. Using the heat kernel techniques presented in Ref. [6] and already 
applied to massless QCD [7], we can obtain the lowest order chiral Lagrangian 



*2 



4 2) = ff-trffarfBpU + tfU + xrf)) (17) 

where U is the non-linear transforming pseudoscalar Goldstone field, x the quark mass 
matrix and try is the trace in flavor space. The pion weak decay constant, /*, at finite 
temperature in the presence of the Polyakov loop and in the chiral limit is given by 

d 3 k 



f*2 = 4M 2 TT V- f 1 

U C kl (27r) 3 [o)2 + k 2 + M 2 ] 



The full calculation of the low energy constants at order 0{p A ) as a function of temperature 
and the Polyakov loop is carried out in Ref. [15]. The main feature is, similarly to (qq)* 
and /*, a strong suppression 0(e~ NcM/3 ) at low temperatures, but an enhancement of 
quark thermal effects close to the chiral-deconfinement phase transition. 

4 Note that in this formula (g(n/3)g(0)) refers to quarks uncoupled to the Polyakov loop while in Eq. 
(8) it refers to quenched QCD. 
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6 Conclusions 



We see that the coupling of the Polyakov loop to chiral quark models at finite tempera- 
ture accounts for large gauge invariance and modifies in a non-trivial way the results for 
physical observables. On the one hand, such a coupling allows to satisfy the requirements 
of chiral perturbation theory at low temperatures, generating a very strong suppression at 
low temperatures of quark loop effects. Nonetheless, the onset of deconfinement through 
a non vanishing value of the Polyakov loop accounts for a chiral phase transition at 
somewhat similar temperatures as in the original studies where the Polyakov loop was 
set to one. We expect this feature to hold also in the calculation of other observables. 
Although these arguments do not justify by themselves the application of these chiral 
quark-Polyakov models to finite temperature calculations, they do show that they do not 
contradict basic expectations of QCD at finite temperature. 

This work is supported in part by funds provided by the Spanish DGI with grant no. 
BMF2002-03218, Junta de Andalucia grant no. FM-225 and EURIDICE grant number 
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